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ANALYSIS OF PRESSURIZED AND AXIALLY LOADED 


ORTHOTROPIC MULTICELL TANKS 

By Robert E. Blum 
Langley Research Center 


SUMMARY 


Expressions for nondimensional stress resultants and displacements are 
derived for multi cell tanks of orthotropic construction loaded by internal 
pressure and axial compression, where the net axial load is assumed to be below 
the buckling value. The present analysis removes an incompatibility appearing 
in a previous solution for isotropic tanks. For isotropic tanks the axial 
stresses obtained from the present analysis differ significantly from previous 
results; however, all the other stresses are in complete agreement. The cor- 
rected axial stresses are presented in the form of curves. 


INTRODUCTION 


Multicell pressure vessels may find application as fuel and oxidizer tanks 
in large launch vehicles. This application represents a departure from the 
conventional circular-cylindrical-shell pressure vessels commonly used in such 
applications at the present time. A sketch of a cross section of an eight-lobe 
multicell tank composed of thin-walled partial circular cylindrical shells and 
radial webs is shown in figure 1. In general, the partial cylinders and radial 
webs may be of unstiffened, stiffened, or sandwich construction. 

A primary motivation for considering the multicell tank is the efficient 
structural use of the radial webs. Whereas the separate slosh baffles required 
in conventional cylindrical tanks constitute largely parasitic weight, the 
radial webs in a multicell tank reduce fuel sloshing while providing additional 
structural strength. Preliminary studies (refs. 1 and 2) indicate that the use 
of multicell tanks in large launch vehicles can provide a significant weight 
saving over the use of conventional cylindrical tanks. This weight saving is 
due not only to the elimination of separate slosh baffles, but also to more 
shallow bulkheads resulting in lighter transition sections between tanks. 

Expressions for the stresses and deformations in a pressurized multicell 
tank of isotropic construction are derived in reference 3 . In that analysis, 
the resulting axial strains in the partial cylinders and radial webs are not 
compatible because the assumption is made that the axial stresses are uniformly 
distributed and equal in the partial cylinders and radial webs. 



The purpose of the present paper is to present an analysis for multi cell 
tanks of orthotropic construction loaded by internal pressure and axial com- 
pression. The multicell tank is assumed to be long enough so that a state of 
uniform strain exists over the cross section; thus, the strain incompatibility 
appearing in reference 3 does not appear in the present analysis. With linear 
cylindrical-shell theory used as the basic mathematical tool, nondimensional 
expressions for stress resultants and displacements are derived, and results 
obtained are compared with results from reference 3 for isotropic tanks. 
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radius of partial cylinders 
coefficients in equation (32) 

thickness of homogeneous partial cylinders and radial webs, 
re spectively 


uniform internal pressure 

displacement in 0 -direction and normal displacement (posi- 
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axial coordinate 


A 1 ; A 2; A 3 

b 1; b 2j b 3 

c l > c 2 

E w,r jEwjX 


5L Cl , 5 l 

pa^ pa^ 
constants in 
constants of 
hoop bending 


Cp ; and -d- e x , respectively 

pa 

equation (32) 

integration in equations ( 22 ) and (23) 
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elastic modulus of radial webs in radial and axial direc- 
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hoop bending-moment, hoop membrane- stress, axial membrane- 
stress, and axial bending-moment resultants, respectively, 
in partial cylinders 

stress resultant in radial webs in radial and axial direc- 
tions , respectively 


•x 


external axial compressive load 


00 

a 

e 0 

€ w,r 
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hoop transverse shear stress resultant in partial cylinders 

half angle between radial webs 
hoop strain in partial cylinders 

radial strain in radial webs 
constant axial strain 
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half arc angle of partial cylinders 
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v 0 


Poisson’s ratios for partial cylinders associated with axial 
stress and hoop stress, respectively 


v w,x> v v,r 


Poisson’s ratios for radial webs associated with axial 
stress and radial stress, respectively 
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Poisson’s ratio for isotropic tanks 
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The multicell tank is assumed long enough so that a state of uniform axial 
strain exists over the cross section. Then, the axial strain e x is constant, 
and the stress resultants acting on elements of the partial cylinders and 
radial webs do not vary in the axial direction. These stress and moment 
resultants are indicated on cylinder and web elements in figure 2. 



Figure 2.- Stress resultants acting on elements of partial cylinders and radial webs. 


Equilibrium Equations 

For constant axial strain, equilibrium equations for the partial cylinders 
are (see, for example, ref. 4): 

Force equilibrium in the 0-direction: 


d0 




0 


Force equilibrium in tbe normal direction: 


( 1 ) 


d0 


+ N 0 - 


pa 


0 


( 2 ) 
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Moment equilibrium: 


- aQ^ = 0 (3) 

Two additional equilibrium conditions are required for the radial stress 
resultants in the webs and the axial stress resultants. Force equilibrium at 
the junction of the partial cylinders and radial webs requires that 





2 N 0 sin (0 


a) + 2 Q 0 cos (e 





Overall force equilibrium in the axial direction gives 



a sin 6 
sin a 


^w,x - P a ^ 


sin (8 - a) sin 0 
sin a 




(5) 


Stress -Resultant -Strain, Strain-Displacement, 
and Reciprocal Equations 

The homogeneous orthotropic stress -resultant -strain equations for the 
partial cylinders are: 


w 0 = K $( € 0 + Vx€x ) 

N x = K x ^x + v 0 € 0 ) 

and the moment -curvature equation is: 

M 0 = " D 0 X 0 

Similar expressions for the radial webs are: 

^w,r = K w,r( e w,r + v w,x e x 


(6) 

(7) 

(8) 
(9) 
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and 
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w,x 


= Ki 


; W,X^€ 


+ V 


w.r^w 




( 10 ) 


When displacements are assumed to he appropriately small, the strain- 
displacement equations for the partial cylinders are (see, for example, ref. 4): 

Hoop strain displacement 


Curvature displacement 


= — 


_ 1/dv 


d0 


- w 




(id 


( 12 ) 


In addition to the preceding equations, the following reciprocal relationships 
must hold among the elastic constants as a result of assuming orthotropic 
construction: 



(13) 


v w,x _ ^w,x 
v w,r Kw^ r 


(14) 


Boundary and Symmetry Conditions 

Three boundary conditions are to he imposed on the displacements in the 
partial cylinders at the junctions with the radial webs: (l) the rotation of 

a normal to the middle surface of the cylinder must be zero, that is, the cyl- 
inder is clamped, (2) the resultant displacement must be in the plane of the 
radial webs, and ( 3 ) the resultant displacement must equal the constant radial 
displacement of the webs. These conditions are expressed mathematically as 
follows : 




= 0 


( 15 ) 


7 



[v cos(0 - cl) - v sin (9 - ajj^g = 


0 


(16) 


a sin 9 
sin a 





( 17 ) 


where , equation 
symmetry of the 


(l6) has been used to simplify equation (17). 
problem is expressed by: 


In addition , the 


w(0) = w(-0) 


(18) 


v(0) = -v( -0) 


( 19 ) 


Solution of Equations 

Equations (l) to ( 12 ) constitute a system of twelve equations in twelve 
unknowns , that is, two displacements, v and wj six stress resultants, 

Nx; M 0 , 00 , Nw,r, and. N w ,xj and four strains, e x , e^, e V}T , and X0. 

Using equation (3) to eliminate Q0 from equations (l) and ( 2 ) and expanding 

the resulting equations for N0 and M0 in terms of v and w by using 

equations (6), (8), (ll), and (12) gives the following sixth-order set of dif- 
ferential equations: 
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d0 2 d0 5 
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where 


( 21 ) 
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Equations (20) and (21) are integrated to give: 


v = C 


( - ^ — + 0 cos j$J + C2 sin 0 


( 22 ) 


w = Ci| 


i + k e 


sin 0,1-k A j . A , _ / 

+ cos 0-0 sin 0 + C 2 cos 0 


1 + k 


+ av x e x (23) 

K 0 


where the conditions of symmetry given by equations (l8) and (19) and the 
clamped condition given by equation (15) have been imposed on the displacements. 


The constants Cq and C2 and the axial strain e x are determined by 
utilizing the boundary conditions (eqs. (l6) and (17)) and the axial equilib- 
rium condition (eq. (5))- In order to accomplish this task, it is convenient 
first to write the expressions for the six stress resultants in terms of the 
solution (eqs. (22) and (23)). The stress resultants can be expressed in terms 
of Cq, C2, and e x by using the strain-displacement equations (eqs. (ll) 
and (12)), the stress-strain equations (eqs. (6) to (lO)), the reciprocal equa- 
tions (eqs. (13) and (l^)), the first differential equation of equilibrium 
(eq. (l)), and the junction equilibrium equation (eq. (4)). When the nondimen- 

sional constants Aq = — — Cq, A2 = ■■■ C2, and A-* = — — e x are introduced, 

pa 2 pa 2 P a 

the nondimen sional displacement and stress resultants are: 


^2 V = Al (rl ~ k ^ S Q n 9 + 0 cos fij + a 2 sin 0 


(24) 
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+ — I — cos 0-0 sin 0 j + A2 cos 0 + v x A^ - 1 (25) 
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(28) 
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(29) 
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2k 
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2 sin(0 - a) 


4k . , K w,x/. 
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L + k 1 K0 y 


1 " V w,rVx) A 3 
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Evaluation of Remaining Constants 

The constants Ap, A 2 , and Aj are determined by using the axial equilib- 
rium condition given by equation (5)j and the two remaining boundary conditions 
given by equations (l 6) and (17). Expanding these three equations in terms of 
Ap, Ag, and Aj gives the following set of algebraic equations: 
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where 
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ap-i = — - — sin 8 + 0 cos 9 - tan(0 - a) ( — - — s ^ n - + — — cos 0-8 sin 0 
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sin 0 cos(0 - a) K w r 
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Solving equation (j2) for A]_, Ag, and successively yields, 
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± jjjg 


+ 2 — sin(0 - a) - 1 
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(35) 
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v x A^ sin(0 


a ) + / — 2 — Bin 9 + 1 - ~k. cos g _ 0 s i D 0 \ 

1 1 + k 0 1 + k J 


sin(0 - a) - ( — - — sin 0+0 cos 0| 
1 + k 


- sin(e - a) 


(1 - - a) * (1 - S’ -ig- sln 6 

(* - |H ♦ *(i - 9 


For a given orthotropic tank, the displacements and stress resultants due 
to internal pressure and axial compression can now he determined from equa- 
tions (24) to (jl) and (33) to (35)- Note that all of the stress resultants 
except the axial stress resultants contain only the constant Ap. Therefore, 
it is convenient to use the expression for Aj (eq. (35) )j t o reformulate the 
axial stress resultants as follows: 


sin 0 sin(0 - a) 


N x / 2 k A \ > ' sin a 

— = V v [l + -== An COS 0] + 

P a x l i + k 1 7 


+ - £ 


E v h sin a 
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“ v w,x|2 sin(6 - a) - 1^— A x sin aj + 




sin 0 sin(0 - a) 


(i - »x)e - ; - rri( Vx ‘ Vw ' x ) Al sin 8 


sin 0 + 2 E x h 
sin a E„ v h.. 


Now, all of the stress resultants can he determined for a given tank hy using 
equations (26) to (29) for N^, M^, Q^, and N w ^ r , and equations (36) and 

(37) for N x and N w ^ x , where Ap is given hy equation (33)* 


A specific case of interest is the homogeneous, isotropic tank. For this 

— 2k 

case it is convenient to write the stress resultants in terms of Ap = Ap, 

1 + k 

h 2 

where k reduces to , to give 

12a 2 


= 1 + Aj_ cos 0 
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-£ = -Ai sin 0 
pa 


( 40 ) 


^ v - ,r = 2 sin(0 - a) - 2 At sin a 
pa* 


(41) 
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sin 0 sin (0 - a) 
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h w sin 9 
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+ 29 


(42) 


= 2v 
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sin (0 - a) - Aq sin a 
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where 


An = 


sin(e -«>-§£ 


sin a + 


h w/ 

CD 

4h\ 

sin 0 


+ COS 0 1 + — — ' 


0 


4k h \sin 0 


+ cos 0 - 


2 sin 0 ’ 
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(44) 


For reasonable values of k = 


12 a 


2’ 


Aq is small and normally may be neglected 


in the expressions for the membrane stresses N 0 , N w ^ r , N x , anc ^ 

A comparison of the present solution with that from reference 3 shows that 
the axial stress resultants and N w x differ; however, the stress result- 
ants N 0 , M 0 , Q 0 , and N w ^ r are in complete agreement. 


N v x . 

The present nondimensional axial stress resultants — and - J — due to 

pa pa 

internal pressure alone are plotted against 0 - a in figures 3j 4, and 5 ^ or 

h w _ . 1 

h ' 2 J 

The terms containing Aq are found to be negligible when computing the stress 
resultants from equations (42) and (43) for these cases. 


an eight-lobe isotropic tank with v = ~, and ^ = 1 , and 2 , respectively. 
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RESULTS AND DISCUSSION 


Expressions for nondimensional stress resultants and displacements are 
given for orthotropic tanks loaded by internal pressure and axial compression. 
Simplified expressions for the nondimensional stress resultants are also given 
for the case of isotropic construction. A comparison of the present results 
for isotropic construction with the results of reference 3 shows that the 
stress resultants , M 0 , Q 0 , and N w ^ r , are in agreement with the previous 

results. However , the axial stress resultants N x and N v ^ x , differ signifi- 
cantly from the results in reference 3 . 


In figures 3* 4, and 5 the present nondimensional axial stress resultants 
N x N w x 

— and due to internal pressure alone are plotted against 0 - a for 


an eight -lobe isotropic tank with 



Y = 1} \> and 2 > 


respectively. 


(It 


should be noted that 0 - a can not be greater than 90 ° due to geometric con- 
straints , i.e., when 0 - a = 90 ° , the partial cylinders are tangent to the 
radial webs at the junctions. When 0 - a, = 0, the multicell tank degenerates 
into a cylinder with radial webs.) For comparison, the axial stress resultants 
for the solution in reference 3> are also presented. Note that the nondimen- 
sional axial stress resultants from reference 3 imply that the axial stresses 
in the partial cylinders and radial webs are equal. When 0 - a = 0 in fig- 
N x 

ure 3 , — from the present analysis is nearly twice the value given by ref- 

pa 


erence 3 ; while 


^w,x 

pa 


from the present analysis is less than half the value 


from reference 3 . The cause of this large discrepancy in the axial stress 
resultants is easily recognized by comparing the expressions for the axial 
stress resultants (eqs. (42) and (43)) with the corresponding expression from 
reference 3 (eq. ( 38 ) in ref. 3)- Eor v = 0 in equations (42) and (43) the 
axial stress resultants agree with the previous results. Therefore, the dis- 
crepancy between the present results and those from reference 3 is due to a 
Poisson* s ratio effect which has been effectively neglected by the failure 
of the previous analysis to enforce compatibility of the axial strains. 


CONCLUDING REMARKS 


Expressions for nondimensional stress resultants and displacements are 
derived for multicell tanks of orthotropic construction loaded by internal 
pressure and axial compression. Simplified expressions for the nondimensional 
stress resultants are also given for the isotropic tank. The present analysis 
removes the axial-strain incompatibility appearing in a previous analysis for 
multicell tanks of isotropic construction. With this incompatibility removed, 
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the axial stresses differ significantly from the previous results; however, all 
the other stresses are in complete agreement with the previous results. 
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